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1.  Introduction 

A  stochastic  process  X  =  (X( t ) ) t€jR  is  self-similar  with  parameter  H  €  IR 

(H-ss)  if  X(c*)  =  c^X(*)  for  all  c>0,  and  has  stationary  increments  (si)  if 
d  d 

X(*+b)-X(b)  =  X(*)-X(0)  for  all  b  €  R,  where  =  means  the  equality  of  all 
finite  dimensional  distributions.  A  real-valued  stochastic  process 
X  =  (X(t))te^  is  symmetric  a-stable  (SaS)  if  all  linear  combinations 
2^  ,a  X(t  )  have  characteristic  functions  of  the  form  exp{-a|0|a)  for  some 
a  >  0.  Here  a  €  (0,2]  and  when  a= 2,  X  is  Gauss ism. 

*\ 

In  this  paper,  we  study  two  disjoint  classes  of  H-self-similar  symmetric 
a-stable  processes  with  stationary  increments  (H-ss  si  SaS  processes).  One 
consists  of  linear  fractional  stable  processes,  which  are  related  to  moving 
average  stable  processes,  and  the  other  consists  of  harmonizable  fractional 
stable  processes,  which  are  connected  to  harmonizable  stationary  stable 
processes. 

In  Section  2  we  give  a  representation  of  a  SaS  process  with  stationary 
increments  in  terms  of  a  stationary  SaS  process  which  is  either  nonanticipating 
or  fully  anticipating.  The  representation  is  shown  for  l<a£2,  but  it  should  in 
fact  be  valid  for  much  larger  classes  than  SaS  processes. 

In  Section  3  we  introduce  the  linear  fractional  stable  processes  Aq  ^(a.b) 
=  (A  u(a,b; t))  whose  corresponding  nonanticipating  (or  fully  anticipating) 

ttlK 

stationary  SaS  process  is  a  moving  average.  We  show  that  when  l<a<2,  to  each 
line  through  the  origin  of  the  parameter  plane  (a.b)  there  corresponds  a 
distinct  linear  fractional  stable  process  (Theorem  3.1).  This  result  implies 
that  the  two  linear  fractional  processes  defined  independently  in  [12]  and  [21] 
are  different. 

In  Section  4  we  introduce  the  complex  harmonizable  fractional  stable 


■„Vi 
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processes  8^  ^(a,b)  =  (0fl  ^(a.b; t))tep.  whose  corresponding  nonanticipating  (or 
fully  anticipating)  stationary  SaS  process  is  harmonizable.  When  l<a£2  we  show 
that  to  each  line  through  the  origin  of  the  parameter  plane  (a,b)  there 
corresponds  a  distinct  harmonizable  fractional  stable  process  (Theorem  4.1),  and 
we  study  their  domain  of  attraction  in  Section  6  (Theorems  6.1  and  6.2). 

There  is  only  one  (distinct)  real  harmonizable  fractional  stable  process, 
namely  Site  0^  ^(1,1),  and  it  is  shown  in  Section  5  that  it  is  not  a  linear 
fractional  stable  process  when  l<a<2  (Theorem  5.1). 

The  only  reason  the  results  in  Sections  3-5  are  stated  only  for  l<a£2  is 
that  the  representation  in  Section  2  is  established  only  for  l<a£2.  If  the 
representation  in  Section  2  is  proved  also  for  0<a£l,  then  the  same  proofs  would 
establish  the  results  in  Sections  3-5  for  all  0<a<2. 

It  should  also  be  mentioned  that,  under  mild  regularity  conditions,  the 
linear  and  the  harmonizable  fractional  stable  processes  introduced  here  seem  to 
be  the  only  self-similar,  symmetric  stable  processes  with  stationary  increments 
whose  corresponding  nonanticipating  (or  fully  anticipating)  stationary  stable 
processes  are  moving  averages  and  harmonizable  respectively.  This 
characterization  result  is  still  under  study. 

The  authors  gratefully  acknowledge  insightful  discussions  on  self-similar 
stable  processes  with  Murad  Taqqu  and  Florin  Avram. 


2.  Representation  of  si  SaS  processes 

In  this  section  we  assume  l<a£2. 

Let  X  =  (X(t))t€jp  be  a  continuous  in  probability  si  SaS  process  with  l<a<2. 
Def ine 

(2.1)  Y(t)  =1°  eu[X( t)-X( t+u)]du  =  X(t)  -  J^e'^'^Xf u)du,  tOR, 

— CO 


A 


where  the  Integrals  exist  a.s.  Clearly  Y  =  (Y { t ) ) is  a  continuous  in 
probability,  stationary  SaS  process  which  depends  linearly  on  the  past 
increments  (or  values)  of  X.  We  call  Y  the  nonanticipating  stationary  SaS 
process  corresponding  to  the  si  SaS  process  X,  and  a  straightforward  calculation 


shows  that  for  all  s<t. 


(2.2) 


X(t)-X(s)  =  Y(t)-Y(s)  -  J-V(v) dv. 


One  can  also  introduce  the  fully  anticipating  stationary  SaS  process  Y 
corresponding  to  the  si  SaS  process  by 

(2.3)  Y(t)  =  j"e~u[X(t)-X(t+u)]du  =  X(t)  -  J‘"e-(v_t)X(v)dv.  t€K. 

and  derive  likewise  representation  (2.2). 

Of  course  the  increments  of  a  stationary  process  Y  define  a  si  process  X 
via  X(t)-X(s)  =  Y(t)-Y(s),  and  so  does  the  indefinite  Integral  of  a  stationary 
process  Y  with  a.s.  locally  integrable  paths  via  X(t)-X(s)  =  J^Y(v)dv.  In  fact, 
a  si  process  is  the  indefinite  integral  of  a  stationary  process  if  and  only  if 
its  paths  are  a.s.  locally  absolutely  continuous;  and  is  both  the  increment  of 
some  stationary  process  as  well  as  the  indefinite  integral  of  some  (other) 
stationary  process  if  and  only  it  its  paths  are  a.s.  locally  absolutely 
continuous  and  its  derivative  process  is  the  derivative  of  a  stationary  process. 
These  two  simple  classes  of  si  processes  are  distinct  (with  nonempty 
intersection)  but  their  union  is  not  broad  enough  to  encompass  all  si  processes 
(as  is  easily  shown  via  examples).  Representation  (2.2)  is  therefore  the  only 
one  generally  available  for  all  si  processes. 

Representation  (2.2)  follows  also  from  Masani’s  representation  of  helixes 
in  Banach  spaces  ([15]),  by  viewing  X  as  a  helix  in  Lp(n,?,P)  for  l£p<a.  It 
should  be  pointed  out  that  (2.2)  defines  a  si  SaS  process  X  for  each  stationary 
SaS  process  Y  even  when  0<a£l,  provided  the  integral  in  (2.2)  is  well  defined; 
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for  necessary  and  sufficient  conditions  see  [6],  Theorem  4.1.  The  conjecture 
here  is  that  (2.2)  would  be  valid  for  0<a£l  as  well,  but  there  may  be  technical 
difficulties  in  its  derivation.  For  the  sake  of  completeness  we  indicate  the 
simple  proofs  of  (2.1)  and  (2.2)  when  l<a£2. 

Proof  of  (2.1).  The  integral  in  the  middle  expression  in  (2.1)  will  exist  a.s. 


eu  E|X(t)-X(t+u)|  du  <  ®. 


(In  fact  this  is  a  necessary  and  sufficient  condition  in  this  SaS  case  ([5]).) 
But  the  stationarity  of  the  Increments  of  X  Implies  that  E|X(0)-X(u)|  grows 
linearly  in  |u|.  The  equality  of  the  two  expressions  a.s.  is  shown  likewise.  □ 


af  (2.21.  From  (2.1)  we  obtain 

•fgY(v)dv  =  X*X(v)dv  -  X^dvj‘^Be‘^v_u^X(u)du 

=  .TgX(v)dv  -  X^duX(u)J^Vue"(v“u)dv 

=  /^(vjdv  -  (e  s-e  t)J®  euX(u)du  -  /t(l-eu  c)X(u)du 
S  s 

=  -(e  s-e  C)X®  euX(u)du  +  X *eU  ^(ujdu 
m  s 

=  Xl,e"(t"u)X(u)du  -  XLe"(S'U)X(u)du 


=  [Y(t)  -  X(t)]  -  [Y(x)  -  X(s)]. 


i.e.  (2.2).  The  interchange  of  the  order  of  integration  is  justified  by 
Fublni’s  theorem  as  for  (2.1).  0 

The  nonanticipating  or  fully  anticipating  stationary  SaS  process  Y 
corresponding  to  the  si  SaS  process  X  via  the  representation  (2.2)  is  continuous 
in  probability.  Every  continuous  in  probability  stationary  SaS  process  Y 
(0<a£2)  has  a  version  of  the  form 
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(2.4)  Y(t)  =  J^a(t.u)dZ(u),  t€R. 

where  Z  has  independent  SaS  increments  and  control  measure  ji,  i.e., 

(2.5)  E  exp(irJfdZ)  =  exp(- |r  |“l|f  |ad|j.) 

for  f  €  La(jx).  (V(t))t^  is  a  strongly  continuous  group  of  isometries  in  La(n) 
and  a(t,»)  =  V(t){a(0,*)}  ([8]).  Since  Y  has  a  measurable  modification  (being 
continuous  in  probability),  the  kernel  a(t,u)  has  a  version  jointly  measurable 
in  (t.u)  ([18]),  and  we  obtain  from  (2.2)  and  (2.4),  provided  all  integrals  are 
well  defined, 

X(t)  -  X(s)  =  /^[aft.u)  -  a(s,u)  -  J*a(v,u)dv]  dZ(u) 

=  •C0[b(t.u)  “  b(s.u)]  dZ(u) , 

where 

b(t,u)  -  b(s.u)  s  a(t.u)  -  a(s,u)  -  I*a(v,u)dv 

=  [V(t)  -  V(s)  -  /V(v)dv]{a(0.*)}(u). 

The  most  important  examples  of  stationary  SaS  processes  are  moving  averages 
and  harmonizable  processes.  Y  is  a  SaS  moving  average  process  (0<a^2)  if 

Y(t)  =  Oi(t-s)dM(s).  ten. 


where  M  =  (M( s) ) Is  a  SaS  motion  (i.e.,  has  stationary,  independent  SaS 
increments  and  Lebesgue  control  measure)  and  h  €  :  =  La(Leb). 

Up  to  this  point  we  have  only  considered  real-valued  SaS  processes.  When 
dealing  with  harmonizable  processes  it  is  natural  and  convenient  to  consider 


complex-valued  SaS  processes,  indeed  radially  SaS  complex  random  variables  will 


suffice.  A  complex  r.v.  XsX^+iXg  is  radially  SaS  if  X^  and  Xg  are  jointly  SaS 
with  radially  symmetric  distribution,  i.e.  with  zsz^+izg,  E  exp{i&ezX}  = 

s 


r 
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E  exptiJZjXj+ZgXg)}  =  exp{-c|z|a}  for  some  c>0.  A  complex-valued  process 
X  =  (X(t))t€^  Is  racially  SaS  If  all  complex  linear  combinations  5^_jZnX(tn) 
are  complex  radially  SaS  r.v.’s.  Y  is  a  harmonizable  stationary  SaS  process 


(0<a$2)  if 


(2.6) 


Y( t)  =  l" eitXdZ(X).  tOR. 


where  Z  has  complex,  independent,  radially  SaS  increments  and  finite  spectral 
(control)  measure  p,  i.e. 


(2.7) 


E  exp(i9lezJ’fdZ)  =  exp(-|z  |aJ*|f  |adp) 


for  all  complex  numbers  z  and  functions  f  in  LQ(p)  ([4]).  (Unless  the  stable 
distribution  of  Z  is  radially  symmetric,  Y  is  not  stationary  ([23]).)  When  p  is 
absolutely  continuous  with  respect  to  Lebesgue  measure,  <f(X)  =  dp(X)/dX  is 
called  the  spectral  density  of  Y.  When  p  is  Lebesgue  measure,  Z  is  called  a 

A* 

complex  SaS  motion  and  is  denoted  by  M.  A  real-valued  SaS  process  Y  is 
harmonizable  if 

Y(t)  =  «e/^ei1:\lZ(\).  tdR. 

where  Z  is  as  above  necessarily  complex.  For  simplicity  here  we  will  consider 
complex  harmonizable  processes,  while  moving  average  processes  will  always  be 
considered  real. 


3.  The  linear  fractional  stable  processes 


In  this  section  we  consider  the  linear  fractional  stable  processes.  Le 
0<H<1,  a.bdR,  and  M  be  a  SaS  motion  defined  in  Section  2.  For  0<aS2,  H^l/a, 
def ine 


(3.1) 


Aa  H(a.b;t)  =  ;_«{a[(t-u)f 1/a  -  (-u)f 1/a] 


+  b[(t-u)^_1/a-(-u)!J"1/a]}  dM(u).  t€R. 


4^ 


with  the  convention  07=0  even  for  tt<0;  and  for  l<a£2,  H  =  1/a. 


(3.2) 
where 

(3.3) 


Aa  jya(a.b;t)  =  aA(t)  +  bM(t),  tQR. 


A(t)  =  (log|t-u|  -  log|u|)  dM(u).  t€lR. 


The  process  A  „(1,0)  was  introduced  in  £21]  and  called  fractional  Lfevy 

motion,  and  the  process  A  „(1,1)  was  introduced  in  [12],  and  called  fractional 

CL ,  n 

stable  process.  A  in  (3.3)  is  the  log-fractional  stable  process  defined  in 
[11].  The  process  A^  jj,  Hj*l/a,  can  be  defined  for  any  a,  0<a£2,  but  the  process 
A(t)  can  be  defined  only  in  the  case  1  <  a  i  2.  It  is  easy  to  see  that  the 
processes  (Aq  ^  (a.b; t) ) t€jR  as  in  (3.1)  and  (3.2)  are  H-ss  si  SaS.  In  the  rest 
of  this  section,  we  assume  l<a<2. 

The  nonanticipating  stationary  SaS  process  Y  u(a.b:»)  corresponding  to  the 

a,  n 

ss  si  SaS  process  A  „(a,b;*)  of  (3.1)  is  determined  via  (2.1)  as  follows:  For 

CL ,  n 


Ya.H(a,b:t)  =  Co  eU  ('Co  *  aC(t^)J  1/0  ~  (t+u~v)+  1/a] 

+b[(t-v)”"1/a  -  (t+u-v)“_1/a]  }  dM(v)  )  du 

=«rZo(CooeU(aC(t-v)+  1/a-(t+u-v)^  1/<Z]  +  b[(t-v)^  1/a-(t+u-v)^  1/a]}du)dM(v) 


(3.4) 


=  X-»ha.b(t_v)dM(v)* 


where 


,  ,  JO  u  ,  r  H-l/a  ,  *H-l/a,  .  r  H-l/a  ,  xH-l/an  , 

ha  b(x)  =  C,  e  (  aCx+  '(x+u)+  3  +  bCx_  ~(x+u)_  ]  ) 


(3.5) 


=  ah^x)  +  bh2(x). 


(3.6) 


,  t  *  H-l/a  -xP  +  v  H-l/a, 

h^(x)  =  x+  -  e  J 0  e  v  dv. 

.  H-l/a  -xr“  -v  H-l/a , 

(3.7)  h^x)  =  x_  ~  e  J x  e  v  dv. 

Hence  Y^  ^(a,b)  is  a  moving  average  of  the  SaS  motion  M,  i.e.  a  linear  process, 
and  we  call  A^  jj(a,b)  a  linear  fractional  stable  process.  ^(a.b)  is  a 
nonanticipating  moving  average  of  M  only  if  h  ,  (x)  =  0  for  x<0,  ,  and  since  for 
x<0,  h^(x)  =  0  but  l^x)  0,  this  is  true  only  when  b=0. 

As  in  (3.4).  the  nonanticipating  stationary  SaS  process  Yq  i/a(a,k) 
corresponding  to  the  (l/a)-ss  si  SaS  process  Afl  ^^(a.b;*)  =  ®A(*)  +  bM(*)  is 
determined  via  (2.1)  as  the  following  moving  average 

(3.8)  Ya  1/a(a.b;t)  =  Oa  b(t-u)dM(u). 
where 

(39)  ga  b(x)  =  ‘•lW  +  ^(x). 

gx(x)  =  log|x|  -  e~xJJ^)evlog|v|dv. 
g^(x)  =  e  x  for  x>0,  =  0  for  x<0. 

Therefore,  Y^  ^^(a.b)  is  a  nonanticipating  moving  average  of  M  only  when  a=0. 

The  linear  fractional  stable  processes  Aq  ^(a.b)  are  indexed  by  the 

parameters  (a.b)  €  1^.  When  a=2,  the  processes  A„  „(a,b)  are  Gaussian,  called 

z,n 

fractional  Brownian  motion,  and  by  comparing  their  covariances  we  see  that  for 
all  a  and  b  with  |a|+|b|  /  0,  the  processes  A0  „(a,b)  are  multiples  of  the  same 
Gaussian  processes  in  distribution.  The  main  purpose  of  this  section  is  to  show 
that  if  l<a<2  each  line  through  the  origin  of  the  parameter  plane  (a.b) 
determines  a  distinct  process  in  distribution  (up  to  scaling  of  course). 

Theorem  3.1.  Let  0<H<1,  l<a<2,  |a|+|b|^0,  and  A  „(a,b;*)  be  given  by  (3.1)  and 

Ci ,  n 

(3.2).  Then  we  have 


Ca>-b>  Aa.H<a*b:‘>  =  Ca>,-b')ia,H(a''b'^' 

where  C  „(a,b)  is  defined  by  (3.10)  and  (3.14)  below,  if  and  only  if 

(i)  a  =  a'  =  0  or 

(ii)  b  =  b '  =  0  or 

(iii)  aa'bb'  /  0  and  a/b  =  a ' /b ' . 

Proof .  Recall  that  for  two  SaS  random  variables  X  and  Y,  by  writing  E  exp(irX) 

d 

=  exp{- |r ja| |X| |“} ,  it  follows  that  X  and  Y  have  the  same  distribution,  X  =  Y. 

if  and  only  if  they  have  the  same  scale  parameters,  ||X||^  =  ||Y||a.  Thus  to 

compare  the  marginal  distributions  of  A  „(a,b;t)  and  A  u(a',b';t)  we  need  to 

d,n  u,n 

compute  their  scale  parameters. 

First  assume  H  /  1/a.  We  find  from  (3.1)  and  (2.5)  that 


|l1a.H(a-b:t,ll“  * 


Itrt  <|a|a*|b|a)  X7||l-v|H-I/“-|v|H-1/“|“dv 

*  jjiaii-»ri/o  -  bivri/ai“dv  > 


(3.10) 


i:  I*' 


For  simplicity  we  drop  the  subscripts  of  a,H.  It  follows  that  for  each  fixed 


C  ^a.bJAJa.bit)  =  C_1(a' ,b‘ )A(a' ,b‘ ; t) . 


where  we  are  assuming  of  course  that  |a|+|b|  /  0  /  la’ |+ |b' | . 

In  view  of  the  one-to-one  relationship  between  the  si  SaS  process  A(a,b;*) 
of  (3.1)  and  its  corresponding  nonant icipatory  stationary  SaS  process  Y(a,b;*) 
of  (3.4), 


(3.11) 


C-1(a,b)A(a,b; •)  =  C-1(a' ,b’ )A(a‘ ,b’ ; •) 


[.*:*»* 


KOI 


i 


is  equivalent  to 


C_1(a,b)Y(a,b;  •)  =  C  ,b’  )Y(a’  ,b  •.  •)  . 


which,  in  view  of  (3.4)  and  Kanter’s  theorem  [9],  is  equivalent  to 

C_1(a,b)ha  b(t)  =  e  C-1(a‘ ,b*)ha.  b.(t-T)  a.e.(t) 

for  some  a  €  {±1}  and  t€IR  (depending  on  a.b.a'.b');  and  this  in  turn  is 
equivalent  to 

C_1(a,b){ah1(t)  +  bhgtt)}  =  e  C-1(a' ,b‘ Xa'hjft-r)  +  b'hgtt-r)}  a.e.  (t) 
in  view  of  (3.5),  which  is  equivalent  to 


C(a‘ ,b' )  ’  C(a, 


Cfa^.b*) 


To  obtain  the  necessity  of  the  last  condition,  take  t  <  min(0,T),  so  that 
h1(t)  =  0  =  h1(t-r)  and  differentiate  C  ^(a.bjbhgft)  =  eC  ^(a' ,b' Jb'hgft-T)  to 
obtain 

bC  (a,b)(-t)  =  e  b  C  (a  ,b  )(T-t)  a.e.  t  <  min(O.T) 

from  which  t  =  0  and  bC  *(a,b)  =  eb'C  *(a’,b')  follow,  and  thus  also 
aC_1(a,b)  =  £a'C_1(a' ,b). 

If  a  =  0.  by  (3.12),  a’=  0  and  (3.12)  is  satisfied  with  e  =  sgn(bb'),  since 
by  (3.10) 

C(0.b)  =  |b|<  ^||l-v|H-1/a  -  |v|H-1/a|a<iv  -  J %  |v|aH_1dv  }. 

Hence  if  a  =  a'  =  0  or  if  b  =  b'  =0,  then  (3.11)  is  satisfied. 


Now  assume  aa’bb’  ?  0.  Then  (3.12)  Is  equivalent  to 


(3.13)  <|=gr. 


Cffa'.b')  &(  a,t 

|a'|a  ’  |b|a 


(^(a'  ,b’ 


-•.v; -.v.vv.v- 

r.  - . . .  ' 


Putting 

i  r00  Mi  iH-1  /a  |  iH-l/aia  , 

A  =  | I 1-v |  -  |v|  |  dv. 

,,  >  rl  II,  iH— 1 /a  |  iH-l/aia  , 

f(x)  =  J0  |  1 1-v |  -  x|v|  |  dv. 

g(x)  =  (l+|x|a)A  +  f(x). 

we  have  from  (3.10),  cP(a.b)  =  |a|ag(b/a)  =  |b|ag(a/b)  and  thus  (3.13)  is 
equivalent  to 

=  *<!>=*&->•  e(|)=s(r-)> 

which  is  equivalent  to  a/b  =  a'/b'  since  g(l/x)  =  g(x)/|x|a.  This  completes  the 
proof  when  H  j*  1/a. 

When  H  =  1/a,  we  find  ||A1/  (a,b:t)||a  =  |t|  cff  .  (a,b)  with 


(3.14) 


where 


c^,  i/a(a,b)  =  2laia  J*i  Uog(v-l)  -  log  v|“  dv 

+  Sq  |a[log(l-v)  -  log  v]  +  b  |“  dv 
=:  |a|a(C  +  f (b/a)), 

f(x)  =  J*Q  |log(l-v)  -  log  v  +  x|“  dv. 


The  rest  of  the  argument  is  similar  to  the  case  H  /  1/a,  using  instead 
expressions  (3.8)  and  (3.9). 


The  linear  fractional  stable  processes  have  the  following  time  domain 
symmetry 


(3.15) 


Aa,H^a,b:*^  =  Aa,H^a,b:_*^ 


In  general,  if  a  self-similar  process  X  =  (X(t))  __  has  a  version  of  the  form 


X(t)  =  J\j;f(t-u)  -  f(-u)]  dM(u) ,  tOR. 


where  M  =  (M(u))u^  has  stationary  and  symmetrically  distributed  increments, 
then  X  has  the  property  (3.15).  Indeed,  by  using  the  si  property  of  M  and  the 
symmetry  of  the  distribution  of  X,  we  have 

X(-)  =  Cj. f(-u)  -  f(-u)]  dM(u) 
d  „ 

=  X^Cf(-v)  -  f( - v)  dM(v) 

d 

=  -X(— )  =  X(-). 

We  now  comment  on  the  sample  path  properties  of  the  linear  fractional 

stable  processes  A  „(a.b;t)  with  0<a<2.  The  fractional  Brownian  motion 
a  •  n 

^2  }j(a’h‘.  c)  ^ias  always  a  samPle  continuous  version.  Kolmogorov’s  moment 
criterion  implies  that  if  X(t)  is  H-ss  si  with  0<H£1  and  E|X(t)|p<“>  for  some 
p>l/H,  then  X(t)  has  a  sample  continuous  version.  Since  A  „(a.b;t)  is  SaS, 

(Z,  n 

E | A  „(a.b;  t)  |P<“>  for  any  0<p<a.  Thus,  if  1/a  <  H  <  1,  then  there  exists  p  such 
CL ,  n 

that  1/H  <  p  <  a.  and  hence  A  „(a,b;t)  has  a  sample  continuous  version;  still 

CL,  n 

the  paths  have  a.s.  nowhere  bounded  variation,  since  H<1  ([24],  Theorem  3.3). 
When  0  <  H  i  1/a,  the  kernels  of  the  stable  integrals  defining  A  „(a,b;t)  have 
singular  points,  which  implies  that  their  sample  paths  are  nowhere  bounded  (see 
[17]  and  [19],  and  for  a  special  case  [13]). 

4.  The  harmonizable  fractional  stable  processes 

In  this  section,  we  introduce  a  new  class  of  complex  ss  si  SaS  processes 

»v 

Let  0  <  a  i  2,  0<H<  1.  a  £  0,  b^O,  a+b  >  0  and  M  be  a  complex  SaS  motion 
(introduced  in  Section  2).  Define 

r»  eltX-l  ,  .1-H-l/a  ,.1-H-l/a, 

(4.1)  0a  H(a,b;t)  =  —  (aA+  +  b\_  )  dM(X).  t€R. 


It  is  easy  to  check  that  0^  ^(a,b;«)  are  H-ss  si  SaS.  When  a=2,  £e  0g  ^(a,b;*) 
is  the  same  process  as  the  fractional  Brownian  motion,  which  is  the  linear 
fractional  stable  process  with  a=2.  This  can  easily  be  checked  by  calculating 
the  covariance  functions  of  those  processes.  However,  when  a  <  2,  9a  ^(a.b;*) 
is  a  new  class  of  H-ss  si  radially  SaS  processes,  as  will  be  shown  in  the  next 
section  for  the  case  1  <  a  <  2.  In  what  follows,  we  assume  1  <  a  i  2  for 
technical  reasons. 

The  nonanticipating  stationary  SaS  process  corresponding  to  0a  ^(a.b;*) 
given  by  (2. 1)  is 

(4.2)  Ya  H(a,b;t)  =  /^eltX  — (aA*"H~1/a  +  b^~H“1/a)  dM(A),  t€R. 

and  is  thus  harmonizable  of  the  form  (2.6)  with 

pA  1  ,  1-H-l/a  A  .  1-H-l/a.  ~  . 

z(A)  =  Jq  (au+  bu-  )  Miu).  AOR. 

In  view  of  the  harmonizabi 1 i ty  of  its  nonanticipating  stationary  process  we  call 

the  processes  0  „  harmonizable  fractional  stable  processes.  It  is  interesting 

(X  •  n 

to  note  that  the  fully  anticipating  stationary  SaS  process  corresponding  to 
0a  H(a,b;*)  given  by  (2.3)  is 


(4.3) 


r°°  itA 


iA  -  1 


(aA*-H  Wa  +  bA*  H_1/a)  dM(A) ,  tOR. 


Both  stationary  processes  in  (4.2)  and  (4.3)  are  (H-l)-ss,  harmonizable,  and 
they  are  identically  distributed. 

The  complex  harmonizable  fractional  stable  processes  0  u(a,b;*)  of  (4.1) 

CL « rl 

o 

are  indexed  by  the  parameters  (a,b)  €  IR+.  We  now  show  that  each  ray  through  the 

2 

origin  of  the  parameter  space  IR+  determines  a  distinct  complex  process  in 

distribution  (up  to  scaling  of  course).  We  also  show  that  the  real  harmonizable 

fractional  stable  processes  Se  0q  ^(a,b;»)  are  all  multiples  of  each  other  in 

distribution,  namely  each  &e  0  u(a,b;*)  is  a  multiple  of  the  process 

cx ,  n 


& 


(44)  *a.H(t):=  *e  0a.H(1*1:t)  *  *e  ' *  ~i\~  |X l  dM(X) ■  teK- 


r«  sin  tXu  il-H-l/a, 


=  -  SZ,  |X I x~"~xx “dMg(X) 


r°°  cos  tX-1  u  1 1-H-l/a , 


where  M^  and  M2  are  the  real  and  imaginary  parts  of  the  complex  SaS  motion  M. 

It  follows  by  [16],  Theorem  7.2,  that  and  Mg  are  real  SaS  motions,  whose 
increments  are  independent  at  distinct  points  X^  ?  Xg:  dM^(X^)  ■H-dMg(Xg),  and 
are  dependent  at  the  same  point  X  in  the  following  manner: 

(dM^X).  dMg(X))  =  2*(dX)1/aR1/2(G1,  G2) 

a/2 

where  R  is  positive  (a/2)-stable  with  E  exp(-rR)  =  exp(-r  ),  r  £  0,  G^  and  Gg 
are  standard  normal,  and  R,  Gj ,  Gg  are  independent.  It  should  be  mentioned  here 
that  0^  jj(l,l;t)  is  also  mentioned  in  the  recent  paper  [20]. 

Theorem  4. 1  Let  l<a£2,  0<H<1,  a^O,  b£0,  a+b>0,  and  0  „(a,b;»)  be  given  by 

Ct,li 

(4.1).  Then  we  have 


(4.5) 


(aa+ba)  1/a0a  H(a,b;.)  =  (a ^a+b,a)'1/a0a  H(a^ '  .b> ;  V ) 


if  and  only  if 

(i) a=a'=0  or 

(ii)  b  =  b'  =0  or 

(iii)  aa'bb'  0  and  a/b  =  a ' /b ' . 
Also  for  all  a,a',b,b'  we  have 


(4.6)  (aa+ba)'1/a«e  0a  H(a,b;*)  =  (a,a+b,a)-1/a*e  0fl  H(a' ,b' : •) . 

Proof.  In  view  of  (2.7)  the  complex  radially  SaS  r.v.’s  /fdZ  and  JgdZ  have  the 
same  distribution  if  and  only  if  their  scale  parameters  ||/fdZ||^  =  I|f|adp  and 
||/gdZ||^  =  /]g|ad^  are  equal.  Thus  to  compare  the  marginal  distributions  of 
0  u(a,b;t)  and  0  „(a'.b’;t)  we  need  to  compute  their  scale  parameter.  For 

CC.rl  CL ,  rl 


,v.v .‘av.v.v.v.v  • 


w.v.v 


r.fti.  Aft  ft  ft.  a  ft. 


simplicity  we  drop  the  subscripts  a,H  in  0  and  Y.  We  have 


. .  .,,a  i  e  -1  ,  N 1-H-l/a  .  L,l-H-l/a, )a 

ll0(a.b;t)l|g  =  J^J  -t—  (aX+  +  bX-  )• 


)|adX 


=  £.1  |  Sin(i|)|a  (Af-")'1  -  bV'1-")-1)  ox 
=  (aa+ba)  2a  /"  |sln(i|)|a  X-(“Htl>  dX 
=  |t|aV'+ba)  j£  I  sin  X|°  X‘(aH  +  ’>  dX 


It  I0"  (a“W)  Da  H. 


Hence  it  follows  that  for  each  fixed  t>0. 


(aa+ba)-1/a0{a,b;t)  =  (a,a+b'“)  1/a0(a' ,b' ; t) 


In  view  of  the  one-to-one  relationship  between  the  si  SaS  process  0(a,b;*) 
of  (4.1)  and  its  corresponding  nonanticipating  stationary  SaS  process  Y(a.b;*) 
of  (4.2).  (4.5)  is  equivalent  to 

(aa+ba)-1/aY(a,b; •)  t  (a,a+b,a)'1/aY(a' ,b’ ; •) 

which,  because  of  (2.7).  is  equivalent  to 

N  N 

(4.8)  ||  2  z  (aa+bar1/aY(a.b;tn)||“  =  ||  2  zn(a'a+b'a)"1/aY(a'  ,b*  ;  tj  |  |“ 

n=l  n=l 


(4.9) 


||  2  z  (aa*baf1/aY(a.b;t  )  1 1‘ 


=  (a“*baf‘  I  2  r.lt"X|“  (a“x;(1'H)‘‘  *  b“x“(1‘H)’1)  (UX2)^2  dX 

n=l 

=  (aa+ba)-1  /*  (aa|  2  z e^Y  +  ba|  2  z  e  ^Y)  C 1  “H } ~ 1  ( 1 +X2 ) "a/2dX 

n=l  n=l 


m 

3 


«o 


rivi 


Hence  if  one  of  conditions  (i),  (ii).  { i i i )  holds  then  (4.8)  holds,  and  so  does 
(4.5).  This  is  the  "if"  part. 


We  next  prove  the  "only  if"  part.  Since  (i)  and  (ii)  are  easy  cases,  we 
only  show  that  (4.8)  implies  (iii).  As  usual,  we  consider  (4.9),  with  N=2. 
Putting  c  =  b/a  and  taking  Zj=l,  Zg=-i ,  t^=0,  t2=t  ±  0,  (4.9)  becomes 

— —  { |l+leitA|“  *  ca|l*le-ltXr>  (l+X2)-“/2  dX. 

l+ca  U 

Therefore  (4.8)  becomes,  with  c'  =  b'/a‘, 

(4.10)  - ic  {|l+ieltX|“  -  |l+ie“*tX|a)  xa(1"H)_1  (l+X2)-®72*^ 

(l+ca)(l-c,a)  ° 

=  2—lc ,a~^l  f*  {(1-sintX)0^2  -  (l+sintX))®72}  Xa( ^^"^l+X2)-®7^  =  0 

(l+ca)(l-c'a)  0 

for  all  t  £  0.  We  now  show  that  (4.10)  is  possible  only  when  c  =  c’.  To  this 
end,  it  is  enough  to  show  that 

f ( t)  :=  X“  {(l+sintX)®72  -  (l-sintX)®72}  xa(1"H)_1  (1+X2)"*172  dX  ^  0 
for  some  t  >  0. 

We  first  consider  the  case  1  i  aH  <  2.  We  have 


f ( t)  =  (Jj/2t  +  j“/2t)  {(l+sintX^-fl-sintX)®72}  Xa{1  H)  1  (l+X2)-®72  dX 
=  J‘j/2t  {(l+ftX)®72  -  (l-ftX)®72}  xa(1_H)-1(l+X2)_a/2  dX 
-  2 j“/2t  Xa(1'H)_1(l+X2)'a/2dX 
>  w  C  S\ D/2t  (l+A2)”"72  dX  -  2  J®/2(.  X-0*1"1  dX. 

In  the  following,  c.(*,»)  will  denote  positive  constants  depending  only  on  the 

J 

parameters  in  the  parentheses.  If  1  <  oil  <  2,  for  sufficiently  small  t>0  we 
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aH 


f(t)  >  Cl(a.H)  t  -  c2(a.H)  t  >  0 


If  aH  =  1.  also  for  sufficiently  small  t  >  0  we  have 


f(t)  >  c3(a)  t  log  -i-  -  c4  t  > 


We  next  consider  the  case  0  <  aH  <  1 .  By  changing  variables,  we  have 

„aH  r°°  , , .  ,  .a/2  .a/2,  a-aH-1  ,  2 ±  2, -a/2 

f(t)  =  t  {(1+sinx)  -  (1-sinx)  }  x  (t  +x  )  dx 

_aH  r®  .  .a/2  , ,  .  .a/2,  -aH-1 

=  t  Jq  {(1+sinx)  -  (1-sinx)  }  x  dx 

aH  r"  .  .a/2  , .  .a/2,  , , .  2  .  2. -a/2  .  -aH-1  , 

+  t  Jq  {(1+sinx)  -(1-sinx)  }  {(1+t  /x  )  -1}  x  dx 

.aH 


(4.12)  =:  t““  (c5(a,H)  +  g(t)>. 


— aH—l 

Since  x  is  strictly  decreasing  on  (0,«).  we  have  c,-(a,H)  >  0.  Choosing 

t  >  0  and  t  >  0  sufficiently  small  we  obtain 

|g(t)|  i  Jq  |(l+sinx)a/2-(l-sinx)a/2|  (1  -  (l+^/x2)-072}  x-0*1-1  dx 


$  2  J*  p®  t  -aH-1 


2e 


0 

1-aH 


"dx  +  2  ;e  ^  x 

X 


dx 


2t‘ 


<  c5(a’H) 


1-aH  (aH+2)e 

Hence  by  (4.12),  we  have  f(t)  >  0  for  sufficiently  small  t  >  0.  This  concludes 
the  proof  of  the  "only  if"  part. 

The  case  of  the  real  processes  Sle  0(a,b;*)  is  easier  to  handle,  because 
their  scale  parameter  is  as  in  (4.7),  and  in  (4.8)  we  need  to  consider  only  real 
coefficients  r  .  in  which  case  (4.9)  simplifies  to 

II  v  f  a  v.a,~l/aai  v  /  ,  .  \  I  la  r“  I  v  ^  *Cn\a  \aO“H)”l  rt  ,%2.-a/2 

II  2  r  (a  +b  )  Se  Y(a.b;t  )||  =  J0  |  I  y  I  X  (1+x  ) 

n=l  n=l 


and  shows  that  (aa+ba)  Y(a,b;*)  =  (a,a+b’a)  ^aSe  Y(a’,b’;*)  for  all 


a.b.a'.b',  and  hence  (4.6). 


0 


By  a  reasoning  analogous  to  that  at  the  end  of  Section  3,  if  1/a  <  H  <  1, 
then  9a  ^(a,b;t)  has  a  sample  continuous  version,  whose  paths  have  nowhere 
bounded  variation.  (When  0  <  H  £  1/a  the  paths  of  0^  ^(a,b;t)  are  expected  to 
be  nowhere  bounded . ) 


5.  The  harmonlzable  versus  the  linear  fractional  stable  processes 

In  this  section  we  show  that,  when  l<a<2,  the  real  harmonlzable  fractional 
stable  process  introduced  in  Section  4  is  different  from  the  linear  fractional 
stable  processes  discussed  in  Section  3. 

Theorem  5 . 1  Let  l<a<2.  Then  the  law  of  the  real  harmonlzable  fractional  stable 
process  ^(*)  of  (4.4)  is  distinct  from  the  laws  of  the  linear  fractional 
stable  processes  A^  jj(a,b;*)  of  (3.1)  when  H  /  1/a,  and  A(»)  of  (3.3)  when 
H=l/a. 

Proof.  We  first  consider  the  case  H  /  1/a.  Recall  from  (4.6)  that  | |#a  ^(t)||a 
=  |t|H21/aD^“  =:  |t|HnraH  and  from  (3.10)  that  |  |Aq  H(a,b;  t)  ||a  =  1 1  \\  H(a.b) 
=  •’  1 1 1  6^  jj(a.b).  For  simplicity  we  delete  the  subscripts  a,H.  We  will  assume 

-1  ^  -1 
ir  X*(-)  =  6  XA(a,b;  • ) 

and  we  will  reach  a  contradiction.  This  implies  the  equality  in  distribution  of 
the  corresponding  nonanticipating  stationary  SaS  processes  given  in  (4.2)  and 
(3.4): 

i  m  j  |\  1 1—  1/a  d  .  _ 

% - <«(X).  t«)  =  (5'V^  b(t-s)dM(s).  t€R). 

Introducing  the  independent  radially  SaS  increments  process  N= 


\r*J**r*  -  A-V'.  "m  m  0  +m+rn.  •  A  • «  •*.  -  .  ..  «*_  /, 


i & 


..  .  |  il-H-l/a  ^ 

N(X)  =  f  Xj  |UL  2  1/2  X  €  R> 

^  °  (i+u 

1  — H— 1  /rr  O  —1/2 

with  spectral  measure  u'-  du(X)/dX  =  ( 6/-r )  | X |  (1+X  )  ,  it  follows  that 

{2*e  C,  eitXdN(X).  t€R}  =  u!»ha  >b(t-s)dM(s).  t€»} . 

Since  the  spectral  measure  N  has  no  atoms,  the  usual  inversion  relationship 
gives  with  I=(x,y),  0^x<y, 


NS(I}  :=  I  {N{I)  +  N{_1)}  =  \  {N(y)  "  N(X)  +  N(_X)  ‘  N(_y)} 

i  t  -iyt  -ixt  .  x 

=  lim  £  J_T  Z—Zfr1 - MO  ^(X)]  dt, 


where  the  limit  is  in  probability  or  in 


-norm.  For  simplicity  we  put 


=  Sje  iutdu  =:  fj(t) 


and  drop  the  subscripts  in  h  ,  .  It  then  follows  that 

SL  •  D 

Ns(I)  t  H™  2^/IT  fj(t)  {  J^,h(t-s)dM(s)  }  dt 
t-t® 

=  lim  27  {  jITfj(t)h(t-s)dt  }  dM(s) . 

Now  ||XfdZ||a  =  ||f||  implies  that  XfRdZ  converges  in  | |*| la-norm  if  and  only 

if  f  converges  in  L  .  It  follows  from  the  above  convergence  in  ||*|l  norm  as 
n  a  & 

T 

T-**>,  that  the  integrand  X_yf  j( t)h( t-s)dt  converges  in  La  as  T-*»,  and  hence  for 
a.e.(s)  along  some  subsequence  Tn  -»  °°.  But 

lim  X^jf  j(t)h(t-s)dt  =  X^f  j ( t )h( t-s )dt  for  a.e.(s) 


by  Lebesgue  dominated  convergence,  since 


f  |  -iyt  -ixt ia'  1/a 

x!jfi(t)h(t-s)|dt  1  X_J- - 1 ir - 1  dt  X_40|h(t-s)  I  dt 
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because  a‘>l  (1/a  +  1/a'  =1).  We  thus  obtain 


1  r“ 


Ng(I)  =  2F  {  ^fI(t)h(t"s)dt  }  dM(s). 


It  follows  likewise  that  the  joint  distribution  of  N  (Ij),  Ng ( Ig)  *s  the  same  as 
that  of 


27  Cn  (  Ol  (t)h(t-s)dt  }  dM(s) ,  <  Oi  (t)h(t-s)dt  }  dM(s) 

1  2 


For  disjoint  intervals  Ij  and  Ig  in  IR+,  the  independence  of  the  former  pair  of 
random  variables  implies  that  of  the  latter.  But  with  l<a<2,  Jg^dZ  and  JggdZ 
are  independent  if  and  only  if  =  ®  a.e.  (in  this  complex  case  see  [4]). 

Therefore  with  1^  =  (0.x).  Ig  =  (a.b),  0<x<a<b,  we  obtain 


_  -ixt  .  m  -ibt  -iat 

- — -7T— -  h(t-s)dt  •  2 - — - h(t-s)dt  =  0 


-it 


a.e. (s) . 


The  same  argument  used  to  show  integrabl 1 i ty  of  fj(»)h(»-s),  establishes  the 
continuity  of  the  integral  f  mf j(t)h(t-s)ds  as  a  function  of  s,  since  the  map 


Has-*  h(*-s)  €  La  is  continuous.  Thus  the  above  equality  holds  for  all  s. 


Putting  s=0,  we  obtain 
-ixt 


(5.1) 


-  1 


-it 


h(t)dt 


«  ^ (b~ a)  _  i  —iat 

- - -77 - -  e  h(t)dt  =  0. 


-it 


But  for  h  €  Lfl  with  Fourier  transform  H  €  Lq.  we  have  (cf.  [22],  Theorem  74) 


.  _  -ixt  , 

H(x)  "  d£  X-»  h(t)  5 dt 


-it 


a.e.  (x) 


and 


,*/  %  d  p®  e  ^  —1  — iat^/^j 

H(y+a)  =  57  -Jt  ~  «  h(t)dt 


(which  follows  by  exactly  the  same  proof  as  for  Theorem  74  in  [22]).  Thus 
differentiating  (5.1)  with  respect  to  x  and  y=b-a  we  obtain  H(x)H(b)  =  0  for  all 
0<x<b,  where  H  is  the  Fourier  transform  of  h.  It  follows  that  H  =  0  a.e.  on  !R+; 


WM 


■%  •*,  .*  v.**,  ■  V*\  •  .  A  «.  •. -‘X.'V.  ■*  A*. -  .A  A  A  A 


and  likewise  on  R_.  Thus  H  =  0  which  implies  h  =  0  a.e.  Because  of  the  linear 
independence  of  h^  and  1^  in  (3.6)  and  (3.7),  h&  ^  =  0  implies  a  =  0  =  b  which 
contradicts  a  +  b  /  0. 

The  case  H  =  1/a  can  be  treated  similarly.  D 

Now  consider  for  l<a£2  the  real  harmonizable  fractional  stable  process 

*a,l/a(t>  =  *e  Xl>  g~4f  1  lx!1_2/a  <«(M.  t€lR. 

which  is  (l/a)-ss  si  SaS  process.  It  is  distinct  from  SaS  motion  (the  simplest 
(l/a)-ss  si  SaS  process),  as  it  does  not  have  independent  increments.  As  we 
have  seen  in  Theorem  4.1  it  is  also  distinct  from  the  log-fractional  stable 
process  A(*),  introduced  in  [11]  as  a  new  (l/a)-ss  in  SaS  process.  Thus 
'J'a  is  a  new  example  of  a  (l/a)-ss  si  SaS  process. 


In  [7],  Theorem  3.1,  it  was  shown  (using  Beurling’s  theorem)  that  the  class 

of  complex  nonanticipating,  invertible,  SaS  moving  average  processes  is  disjoint 

from  that  of  the  regular,  harmonizable,  stationary  SaS  processes.  The  proof  of 

Theorem  5.1  shows  that  the  two  classes  of  processes,  {^eJ‘_#oe  dZ(X),  t  €  !R}  and 
00 

{/  h( t-s)dM(s) ,  t  €  R},  are  disjoint.  Namely,  the  entire  class  of  real 


harmonizable  stationar 


average  processes. 


rocesses  is  disjoint  from  that  of  real  SaS  movi 


6.  The  domain  of  attraction  of  the  harmonic  fractional  stable  process 


The  domain  of  attraction  of  the  linear  fractional  and  of  the  log-fractional 
stable  processes  has  been  studied  in  [1],  [2],  [10],  [11],  [12].  In  this 
section,  we  study  the  domain  of  attraction  of  the  H-self-simi lar,  symmetric 
a-stable  process  with  stationary  increments  0  „(a,b;*)  introduced  in  Section  4. 

CL ,  n 

We  start  with  0  „(1.1;*). 

a ,  n 
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Theorem  6.1.  Let  1  <  a  <  2,  1-1/a  <  H  <  1,  M  be  the  complex  SaS  motion  in 
00 

(4.1),  (Yn)n _ m  be  the  harmonizable  SaS  sequence 

Yn  =  einX  dM(X) 


(6.1) 


Z  =  2  c  ■ Y. 

n  i  m  n-k  k 
k=-°° 


with  appropriate  coefficients  (c  }.  Let  t  =  H+l/a-1  €  (0,1).  Define  {c  }  by  c„ 

n  v  n  u 

=  0.  cn  =  |n|T  *  for  n  /  0.  Then  (Zr}  is  well-defined  by  (6.1)  and  as  n-f“>, 

i  Cnt]  d 

“H  2  Zm  CTfla.H^1,1:t^*  *  >  °* 

n  m=l 


where  =  4T(tt)  cos(-nr/2),  and  -»  means  the  convergence  of  all  finite 
dimensional  distributions  of  the  Indicated  process. 


Proof ■  (Step  1)  We  show  that  (Zn)  is  well-defined.  Since 


?  cn_kYk  =  einX  (  2  cm  e'lmX)  dM(X). 


|k|SK 


m=n-K 


it  suffices  to  show  that  the  Fourier  series  2°°  c  e  (i.e.  the  sequence  of 

m=-®  m  v 

its  partial  sums)  converges  in  L^-w.ir).  Since  -r  <  1,  c  =  |m|  1  0  as 

|m|  t  00  and  thus  2^__<ncm  e  is  convergent  everywhere  except  at  X  =  0  (mod  2 it) 
(see  [3],  pp.  87-88).  Also  since  a  >  1,  a  necessary  and  sufficient  condition 
for  2  c  e  €  L  (-w.ir)  is  2”  _ca  |m|a  2  <  ®  (see  [3],  p.  207),  which  is 
satisfied  since  a(t-l)+a-2  =  a(H-l)-l  <  -1.  In  fact  the  argument  in  [3],  p. 

208,  shows  that  in  this  case  the  sequence  of  partial  sums  converges  in  L^f-ir.ir). 
For  completeness,  we  show  this  fact  below. 

For  simplicity,  we  consider  the  convergence  of 

m  00 

f  (x)  :=  I  c,  cos  kx  to  f(x)  :=  2  c,  cos  kx 

m  k=l  K  k=l 


in  L  (O.ir)  as  m  -»  «.  By  the  same  argument  as  in  [3],  p.  208,  for  n  >  m,  if 
7r/(n-nH-l )  <  x  <  ir/  ( n-m ) , 

n  » 

|f  (x)  -  f (x)  |  12  c.  +  |  2  c,  cos  kx| 
k=nH-l  k=n+l 

n  n  n 

12  c,  +  ^  c  1  2  c,  +  const. (n-m)c  12  2  c,. 

k=m+l  k=m+l  k=m+l 

Hence,  with  Bn  =  ck. 

J'S  |fm(x)-f(x)l“dx  -  If„(x)-f(x)l“-ix 

<  const.  2  ^/(n-ie+l)  Bn  5  const'  JL,  Bn  <n-|”>~2 

n=tn+l  *•  ’  n=m+l 


o 

1  const.  2  B®  i 
«=1  *+m 


by  a  similar  argument  as  in  Theorem  2  in  Appendix  22  of  [3], 

00  00 

1  o  (1)  +  const.  2  c“  €a~2  1  const.  2  for  large  m 

m  £-1  ^+m  «=1 

We  will  show  lim  _  l“  ,  c“  €a~2  =  0  under  our  assumptions.  Let 
m-*»  g=l  e+m 

A  :=  2  c“  e°~2  <  “. 

«=1  E 

For  any  e  >  0.  there  exists  L  =  L(e)  such  that  c^ia  2  <  e.  For  any 

.  .  .1/(1-- r)a, 
m  >  LA  v  ’  /t, 

00  co 

-r  “  oa_2  ✓  -e  «  »«-2  w 

2  C-  €  ^  2  C«  C  <  t. 

«=L+1  *+m  £=L+1  * 


Note  that 


-r-  i  -  <  -T77T — T“  for  e  *  L- 

£+m  m  ^1/(1-Tr)a 


Hence 


-by 


1—7  1-T 


Therefore  any  m  >  LA 


l/(l-i)a/ti 


_  a  na-2 


(l-'r)a  L 


2  ca  ea~2 
e=i  e 


(l-Tr)a 


Hence 


.  .  _  a  .a- 2 

lim  sup  I  c,,  6 

^  n  -  e+m 

mH»  5=1 


£  +  £ 


(l-Tf)a 


and  so 


lim 


00 

2 


m-*»  e=i 


ca  ea~2  =  o 
ce+ m  *  u* 


concluding 


lim  S0  _  f^x^!a  dx  =  0. 

nn® 


(Step  2)  We  have  for  each  t  >  0, 


Z  (t) 
nv  ’ 


n 


[nt] 

2  Z_ 
m=l 


m 


[nt] 


1  rir  ,  imX>  ,  v  -ikXA  ,1.  > 

=  ji  J_T  (  2  e  )  (  2  e  )  dM(X) 


n 


1 


m=l 


r«mr 


k=-« 


imu/n 


nr-l/a+1  -rnr 


i  [nt] 
rnir  ,1  L_  J 
=  J  (—It 
-rnr  vn  . 
m=l 


(  2  e 

m=l 

imu/n 


)  (  2  c  e 
k=-°° 


-iku/n 


)  dM(^) 


)  (^  2  c^e”^ku/n)  dM(u) 
n  k=-°° 


(6.2) 

where 


=  ^nw  Kn(t’u)Lr>>  =:  Zn(t). 


n 


n' 


v  , .  ,  1  ^  imu/n 

K  (t,u)  =  —  2  e 

nv  ’  n  , 
m=l 


1  iu/n  elu[nt]/n  _  i 


=  —  e 
n 


n-n® 


e 

itu 


iu/n 


-  1 


iu 


|u  |  i  nir, 

0  <  |  u  |  £  nir , 

u  /  0, 


>5 
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S-lir  lKn<t*u)|a  lLn(u)  "  Ln(u)  1“  du  -  °- 


For  u  /  0  we  have 


|Ln(U)-Ln(u,|  =  4  |  I  (c,  -  =i).-‘Wl’l  1 


n  k=-°° 


since  by  (6.3),  c^  -  c^  =  0(|k|  ^).  Also  for  0  <  |u|  £  ror, 


(6.10)  |K  (t,u)  |  =  — 


sin(j£) 


2  sln(IntJu}  _2n. 


Sin%) 


Since  nr  >  0  and  a  >  1,  it  follows  that  with  0  <  a  <  nir. 


f  |K  (t,u)|a  |L  (u)  -  L‘ (u) |a  du  < 


const. 


a<  |u  |<mr 


nay  a<|u|<mr  |u|a 


(6.11) 


Thus  it  remains  to  show  that 


const. 


0  as  n  -*  °°. 


X  |K  (t,u)|“  |L  (u)  -  L'(u)|a  du  -w  0  as  n  -»  ». 

|u|<a  n 

However,  by  the  argument  in  [14]  (see  Equations  (2.10),  (2.16),  and  the  next  one 


on  p.  139  in  [14]),  we  see  that 


1  |u|<iw  IV'-")'2  IL„<U>  -  Ln<“>|2  lulq''  du  -  °  "  " 


for  any  q  €  (0,1).  Hence 


/ |u |<a  lyt.u)!2  |Ln(u)  -  Ln(u)|2  |u|q_1du  0  as  n 


which  implies  (6.11)  by  Holder’s  inequality. 


(Step  41  We  finally  show  that  for  each  t  >  0, 


(6.12) 


Zn(t)  ->  C70H(l.l;t) 


5 


|| 
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in  probability,  or  equivalently  that 


(6.13) 


~  itu  - 

llz^o^e  -  C,  I  Vt  u)Ln(“)  -  1“!“'  1°  d“ 

itu  . 

+  S  |  C  6  — — |u  |  |a  du 

|u|>ror  1  lu 

-*  0  as  n  ■*  ®. 


The  integral  J*|u|>nir  tends  to  0  as  n  ■*  “  because  the  integrand  in  Lebesgue 
integrable  over  R.  For  the  integral  over  |u|<mr  we  have  from  (6.6), 

itu  _  . 

gn(t.u)  :=  KJt.u)  1/  (u)  -  ■  |up 

itu  . 

=  C  lup  {  K  (t.u)  f(-)  -  - — }. 
nr 1  1  1  nv  ’  ny  in  1 

It  then  follows  from  (6.7)  that  as  n  ■*  •, 


gn(t.u)  -*  0.  u  *  0, 


and  from  (6.8)  and  (6.10)  that 

|gn(t.u)|  i  Cju|  ‘T  min  (1+t,  jjjj-  },  0  <  |u|  i  mr, 

where  the  function  on  the  right  hand  side  Is  in  La(R)  since  err  =  l-a(l-H)  <  1 
and  a(-r+l)  =  aH+1  >  1.  Hence  the  dominated  convergence  theorem  implies 

J"  |g  (t.u)  |a  du  -»  0  as  n  -»  <® 

euid  thus  (6.13)  is  established. 

It  then  follows  from  (6.12)  that  as  n  ®, 


Ct  Jf1aj0a.H(1,1:tj) 


in  probability,  and  the  result  follows  from 


m 


bw 


which  is  established  as  (6.2),  and  from  (6.9). 


The  restriction  1-1/a  <  H  <  1  in  Theorem  6.1  was  used  only  in  the  estimate 
(6.11).  If  the  left  hand  side  of  (6.11)  could  be  shown  to  converge  to  zero  for 
0  <  H  <  1,  then  Theorem  6.1  would  hold  for  all  0  <  H  <  1. 

With  some  modification  of  the  proof  of  Theorem  6.1  along  the  idea  of  Major 
[14],  we  obtain  the  following  generalization. 

Theorem  6.2  Let  a£0,  b£0,  a+b>0  and  define 

A  =  -  I  a+b  -  i  a~b  i 

8T(t)  '  cos(ttt/2)  sin(nr/2)  '* 

where  t  =  H+l/a-1  as  In  Theorem  6.1.  In  the  assumption  of  Theorem  6.1,  we 
replace  c^  =  In^  * ,  n  /  0.  by 


An7  * .  n  >  0, 

n  1  A|n |”r_1 ,  n  <  0. 


Then  we  have 


1  tnt-l  d 

7  0«.H(a-b:t>- 

Sketch  of  Proof.  We  give  only  the  out  ling  of  the  proof  here  and  omit  the 
details. 

With  the  same  function  f.  as  in  the  paragraph  preceding  (6.3),  we  define. 


instead  of  the  {c^}  of  (6.3), 


^  :=  §7  K^jreikX|Ar7f(X)dA.  \  K;J^e1KA|\pf(X)  sgn  X  dX 


•  .  i  •  i*w  _ i  kX  | .  |  — r . 


where 


=  2(A+A)r('r)cos(nr/2) ,  IT  =  2(A-A)r(-r)sin(-nr/2) 


WS5SRK 


f ft.J 
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Then  by  the  same  argument  as  in  [14], 


a^  =  ±(A+A)  Ikf'1  +  0(k~*).  b^  =  i(A-A)|kr  'sgn  k  +  0(k"*) 


-2 


1  /  i  A  \  It.  I  1 


-2, 


‘r_1  +  0(k-2) ,  k  >  0, 

Alkf-1  +  0(k-2) .  k  <  0, 


implying 

f  A  |k  | 

(6.14) 

l< 

n 

jf* 

+ 

rf* 

n 

d* 

and  also 

2  aJ  e  ikx  =  K  |x|  ^f(x).  2  b/  e  ikx  =  iK' |x|  ^f(x)  sgn  x. 

k=-«  k=-»  K  y 


Now  put 


with 


2n(t)  *  ^nw  Kn(£*u)Ln(u)  dM(u). 


.  -iku/n 


L„<“>  -  h  ,L>  • 

n  k=-« 


=:  LA.l(u)  *  l;.2(u>- 


0  <  u  £  rnr. 


where 


(6.15)  Ln  j(U)  =  Kju|  Tf(J).  Ln  2(u)  =  iK;|u|  sgn  u. 

_2 

Note  that  from  (6.14),  c^-c^  =  0(k  ).  So  the  same  argument  as  in  Step  2  of  the 

proof  of  Theorem  6.1  concludes  that  as  n  -»  *, 

X-nir  lKn(t  u>  1“  lLn(u>  “  LR(U)  T  du  -  0. 

Step  4  will  be  handled  as  follows.  First  note  that 

K^^a+b),  K;  =  -|(a-b), 

so  that 


£y< 
y.y. 

sfci 

ifi 


XA 


Hence  we  have 


a  =  K  +  iK' .  b  =  K  -  IK1 

7  7  7  7. 


m  itu  , 

j*>  e  -1 


0a.H(a>b:t)  =  C.  HT-  ((ViK;>ur  +  (K,-iK;)u_,>  dM(u) 


_  itu  . 

jo  e  -1 


{K>P  +  iK;|upsgn  u}  dM(u) . 


_  __  itu  . 

ll2^*)  '  0o.H^a,b;t^U  =  ^-mr  I  Kn(c-u)Ln(u)  ~  "  iu  * lu l~^(K.Y+1K^sgn  u)  1°  du 


itu  . 

+  x iu|>nw 1  rur~|u|  'Y(K7+lK;s«n  u)  1°  du- 

The  second  integral  above  tends  to  zero.  As  to  the  first  integral  we  have  from 
(6.15)  that  the  integrand  without  the  a  power  is  equal  to 

itu,  itu  . 

Kn(t-u)Ln.i<u>  ■  Krnr"|u|_nr  +  Kn(tiU>LA,2(u) "  iKiLirlu|  Tsgn  u 

itu,  itu 

-  K,|uf  {  Kn(t,u)f(^)  -  }  +  iK^upsgn  u  {  Kn(t.u)f(£)  -  } 

By  the  same  reasoning  as  in  the  proof  of  Theorem  6.1,  the  right  hand  side  above 
tends  to  zero,  and  the  integral  itself  also  converges  to  zero.  □ 
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